Thermal effects on atomic friction 
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We model friction acting on the tip of an atomic force microscope as it is dragged across a surface 
at non-zero temperatures. We find that stick-slip motion occurs and that the average frictional 
force follows | lnu| 2 / 3 , where v is the tip velocity. This compares well to recent experimental work, 
permitting the quantitative extraction of all microscopic parameters. We calculate the scaled form 
of the average frictional force's dependence on both temperature and tip speed as well as the form 
of the friction- force distribution function. 
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The study of friction between two surfaces presents 
numerous theoretical and experimental challenges 
Two macroscopic surfaces interact through many asper- 
ities. These asperities are on an atomic length scale, 
so macroscopic friction is inextricably linked with micro- 
scopic properties. A friction force microscope, the tip 
of which is in essence a single asperity Hf||, provides 
an ideal way to probe atomic friction by removing the 
complication of multiple asperities. Under a constant 
load, the low-velocity motion of the tip on a surface ex- 
hibits stick-slip behavior as the moving tip hops over an 
atomically-defmed substrate pj|-|ic|]. In principle, such re- 
sults can give a quantitative description of friction at the 
microscopic level. However this requires a theoretical un- 
derstanding of the behaviour of the tip on the substrate. 
This is the purpose of this Letter. 

Recent studies have argued that, due to thermal fluc- 
tuations, the lateral friction force F has a logarithmic 
dependence on velocity F cx const + lni> where v 
is the velocity of the support. This corresponds to lin- 
ear creep between two surfaces in contact when the force 
acting on a substrate produces a small constant poten- 
tial bias jllj . But this is quite different from what occurs 
when a tip is dragged across the substrate: the potential 
bias is continuously ramped up as the support is moved. 
We will call this ramped creep and will show that, at 
constant temperature T, ||l^-|l7|| 

Fcx const -T 2/3 | In v/T| 2/3 . (1) 

It is straightforward to discriminate between linear and 
ramped creep. Consider Figs. [j]a and ^b where the com- 
bined potential of the atomic substrate and the elastic tip 
is shown. Locally, it is sufficient to consider two wells. 
In that useful analogy to (mean-field) nucleation 

at a first-order phase transition occurs (l^Jf^l ■ When the 
bias is infinitesimally small, the energy, up to an additive 
constant, is well described locally by E = Fx + x 2 — x 4 , 
where F is friction and x is displacement. Hence the bar- 
rier height is AE ~ (const — F), as for mean-field nucle- 
ation near the coexistence curve. In the steady state, the 
rate of thermal fluctuations is proportional to the veloc- 
ity. Since the rate is itself proportional to the probability 



of a fluctuation over this barrier, one obtains linear creep: 
the force depends logarithmically on the velocity. This 
is however an unlikely regime because the barrier is so 
high (see Fig. 0a). In contrast, as the tip is dragged over 
the substrate and the bias is ramped up, barrier-hopping 
fluctuations occur preferentially when the tip is close to 
slipping at the top of the barrier. At this point, the local 
combined potential can be written E = Fx — x 3 , where 
the origin of the x axis has been conveniently shifted, 
and the second well has been replaced by an absorb- 
ing boundary condition (see Fig. 0b). Then we have 
AE ~ (const — i^ 1 ) 3 / 2 and Eq. (1) for ramped creep is 
obtained. This is similar to mean-field nucleation near 
the spinodal |O,0|, and is calculated, as well as observed 
experimentally, for, e.g., magnetic- flux fluctuations in su- 
perconducting quantum interference devices [jf4|-[f6| , and 
in other contexts |l7| ]. We give a detailed exposition of 
the algebra below. First, however, we solve a model nu- 
merically and show that the experimental regime corre- 
sponds to our physical picture. 

The system composed of the tip of the force micro- 
scope interacting with a surface can be described by a 
modified Tomlinson model {j]-[l(J. We thus consider a 
tip with longitudinal and transverse coordinate positions 
x and z and effective mass M. The tip moves under 
the influence of the surface adiabatic potential E(x,z), 
coupled elastically to the support of the microscope, of 
coordinate R(t) and Z(t). For constant load, the distance 
Z(t) — z is fixed, so the transverse coordinate plays no 
role, and R = vt with v the constant speed at which the 
tip is dragged. Including the random effect of substrate 
fluctuations, we have the Langevin equation 
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where M is the mass of the tip, 7 is the microscopic 
friction coefficient, E is the combined surface-tip poten- 
tial, and £ is a random noise satisfying the fluctuation- 
dissipation relation (£(*)£(*')) = 2M~fk B T5(t—t'), where 
the angular brackets denote an average, and ks is Boltz- 
mann's constant. The combined surface-tip potential has 
the form 



f 



E{R,x) 



£ (R(t) - x) 2 -U cos (— 
2 \ a 



(3) 



where U is the surface barrier potential height, a is the 
lattice constant, and k is the elastic constant of the tip- 
support coupling. There are a a series of potential wells 
x(R) given by the solution to dE(R,x)/dx\% — 0, with 
the instantaneous lateral friction force on the support 
F(t) — k (R(t) — x(t)). A typical sequence of the stick- 
slip motion of F is shown in Fig. lc. 

The relevant time scales of the problem are 7 -1 , the 
viscous time scale for dissipation of energy from the tip 
to the surface [M, po = y/ Ma 2 /U, the oscillation pe- 
riod of the mass in the surface potential, and the tip 
resonance period pu = y/M/k. These parameters are ef- 
fective in the sense that their precise value may depend 
on the load applied to the friction force microscope [jl9| . 
The friction coefficient must be larger than the critical 
value needed to observe overdamped stick-slip behaviour 
at zero temperature. Within the Tomlinson model, this 
value is ~ 2uj k [Q||[l0| ■ To first order in the deviations of 
the tip from an integer position, the instantaneous force 
can be approximated as F(t) ~ [k/(l + f2| )]iJ(t) = kR, 
where tt k = po/C2np k ). 

Equation (||) was simulated using Ermak's algorithm 
pof with parameters consistent with experimental con- 
ditions. The lattice parameter a — 0.4 nm and the ef- 
fective spring constant k = 0.86 N/m [ElJ were taken 
from the data of Ref. || , obtained under constant load 
F n = 0.65 nN. The parameters U = 0.27 eV and 
7 = 8.9 x 10 5 sec" 1 and spring resonance frequency 
(2npk)~ 1 = 52MHz were adjusted to obtain the best 
fit with the experimental data. The mass M — 8.7 x 
10~ 12 kg then follows. The tip oscillation frequency 
2-kPq 1 = 1.1 x 10 6 sec _1 is a comparatively small at- 
tempt frequency compared to usual values in adatoms' 
surface diffusion J22|, which is why thermally activated 
events with very small barrier height are important. 

The numerical simulations were performed over a large 
range of velocities 5 nm/sec < v < 256 /im/sec, with sev- 
eral temperatures 53 K < T < 373 K and time steps 
0.001 po < St < 0.01 po- From the stick-slip motion of 
F, we extracted the average lateral force as a function of 
velocity, for different temperatures. This is shown in Fig. 
[j]d. To test the prediction for ramped creep, we first ex- 
tracted the constant force in Eq. ([!]), F c , from the data 
at different temperatures, but fixed ratio of the scaled 
velocity v/T, as shown in the top inset of Fig. 2. With 
this constant removed, the prediction from Eq. (1) cor- 
responds to a universal form, which is independent of 
temperature. Fig. 2 shows a very good scaling collapse 
of the dependence of the average force on velocity for 
numerical data from five different temperatures, as well 
as for experimental data from the work of Gnecco et al. 
pj . This collapse well confirms the prediction of ramped 
creep from Eq. (1). The small deviations from the loga- 



rithmic behaviour at large values oiv/T> 1 indicate the 
approach to the zero-temperature stick-slip motion 
The bottom inset of Fig. || shows that, in contrast, linear 
creep does not give a good scaling collapse. 

These results can be explained analytically by consid- 
ering the thermal activation transition of the tip out 
a single metastable well and the distribution of the 
support's position P(R) when such a transition occurs 
Jl4] , ^6[ . Translational invariance allows us to concentrate 
on one such minimum, of value x = if R — 0. As 
the support is moved, the barrier potential to the next 



minimum vanishes at some critical position R c 14 1, de- 
termined by dE/dx c — d 2 E/dx 2 c — 0: 
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where 2irx c /a = cos _1 (— fi 2 ). Note that the support's 
motion corresponds to a continuous linear ramping of 
the potential. Provided that the velocity of the support 
is slow enough, there can occur thermally activated tran- 
sitions between two minima before the critical position 
i? c is reached. This is described by the Kramer's rate 
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where f2 and AE correspond respectively to the instan- 
taneous effective oscillation frequency and barrier height. 
Thermally activated transitions are most likely to occur 
when the support's position is close to the critical value. 
To lowest order in the bias /(*) = 1 - (R(t)/R c ) < 1, 
the minimum and maximum of the potential are easily 
obtained, and one finds MM: 
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Note the barrier height vanishes as / 3 / 2 as anticipated 
above. Likewise, the factor f 1 ^ 4 justifies the high-friction 
limit (7 ^> f2) of the prefactor in Eq. (||). 

The transition rate represents a measure of the time 
needed before a thermally activated transition takes 
place, in particular, W(R(t)) — exp — f T^ 1 (R(t')dt' 
gives the probability that a transition has not taken 
place at time t. Due to the exponential character of 
the transition rate, the probability for transition is much 
greater when / <C 1 and insensitive to the initial support 
position. The distribution of the support's position at 
which a transition occurs is then conveniently expressed 
in terms of the reduced bias /* = (AE/k B T) 2 ^ as pi : 
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where the dimensionless velocity 
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is essentially a function of the temperature, microscopic 
friction coefficient and velocity of the support. The av- 
erage of the distribution (/*) is then found from Eq. (||), 
with value (/*) = |lnu*| 2 / 3 + 0(l/ln«*) in the limit 

v* « i 0g. 

Once the average support's position at which a slip 
event occurs is known, it is straightforward to calculate 
the average lateral force as the integral of the instanta- 
neous force over a cycle of the stick-slip motion (24|] ; the 
form of Eq. (1) then follows: 



F = F C - AF\ lnu*| 2/3 , 
where the constants pal 
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With the parameters used in Fig. 2, Eq. (|9|), Eq. ( |ll| ) 
and Eq. Q predict values of F c = 0.55 nN, AF/T 2 / 3 = 
1.91 x 10~ 3 nN/K 2/3 and Tv*/v = 3.9 x 10~ 2 K-sec/nm. 
This compares well to the respective values of (0.54 ± 
0.01) nN, (1.82 ±0.02) x 10~ 3 nN/K 2/3 and (4.5 ±0.2) x 
10 _2 K-sec/nm extracted from the numerical and exper- 
imental work, confirming that the ramped creep regime 
well describes the motion of the tip. 

From the distribution function, other properties can be 
calculated. In particular, a directly accessible experimen- 
tal quantity is the fluctuation of the maximal force F m , 
the force just before a slip event occurs. Rewriting the 
bias as / = l — F m /(kR c ) allows us to calculate the distri- 
bution of the force maxima. The analytic form is shown 
in Fig. 3, together with numerical data. The agreement 
between the numerical results and theoretical prediction 
is very good. Finally, we note that it is also straight- 
forward to redo the above calculation for the restricted 
regime in which linear creep takes place. This regime 
is limited to very small velocities, less than O.lnm/sec: 
an expansion of E(R, x) around x — and x = a yields 
the condition v < (4_ira/~jpl)(k B T /k) exp(-f c U/k B T), 
where f c > 0.5 is the value of the bias at which a linear 
expansion becomes inappropriate. 

To test our results experimentally, it would be partic- 
ularly valuable to consider an extended range of temper- 
atures and velocities, obtaining not only the dependence 



of the friction force on those quantities, but the distri- 
bution function of such forces as well. The parameters 
extracted from such comparisons would provide direct, 
feasible, and simple access to the fundamental descrip- 
tion of friction at the atomic scale. 
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FIG. 1. (a) Schematic representation of the combined 
tip-surface potential just after a slip event. At this point, 
the energy barrier is high and thermally activated transitions 
improbable, (b) As the support is moved, the energy bar- 
rier diminishes and transitions become more likely. Inset (c) 
shows typical stick-slip behaviour of the instantaneous friction 
force as the support is moved at velocity v = 25nm/sec for 
a temperature T = 293 K. Inset (d) shows the average fric- 
tion force for different temperatures. Open circles, squares, 
diamonds, triangle-ups and triangle-downs correspond respec- 
tively to temperatures T — 53, 133, 213, 293 and 373 K, closed 
circles are the experimental data of Ref. [||. The units of 
velocity v are nm/sec. 

FIG. 2. Justification of the scaling form, Eq. (|l|). The 
upper inset shows the value F c as extracted from the data of 
the friction force shown in Fig. [j] for different temperatures 
between 53K and 373K and fixed ratio T/v = lK/(nm/sec). 
The lower inset shows significantly worse scaling for linear 
creep, using ln(v/T). The units of velocity v are nm/sec tem- 
perature is in Kelvin. 



FIG. 3. Normalised distribution, with statistical error 
bar, of the maxima of the friction force F m at temperature 
T = 293 K and velocity v — 25 nm/sec. The result from the 
simulation is well reproduced by the theoretical distribution, 
Eq. (§). 
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